Product year selection

ir_pop <- read_csv(str_c(datadir, "processed/ir_population_data.csv"), show_col_types = F) 

selected_products <- NULL
selected_years <- NULL

# Get product year combinations
plot_specs <- 
  ir_pop %>% 
  select(product, year) %>% 
  unique() %>% 
  mutate(ir_pop = list(ir_pop), .before = everything()) %>% 
  {if (!is.null(selected_products)) filter(., source %in% selected_products) else .} %>% 
  {if (!is.null(selected_years)) filter(., year %in% selected_years) else .}

Density plots

Difference of UN Country Population and Product Country Population

## [[1]]

## 
## [[2]]

## 
## [[3]]

Ratio of UN Country Population and Product Country Population

## [[1]]

## 
## [[2]]

## 
## [[3]]

Maps

## Spherical geometry (s2) switched off
## Warning in st_simplify.sfc(st_geometry(x), preserveTopology, dTolerance):
## st_simplify does not correctly simplify longitude/latitude data, dTolerance
## needs to be in decimal degrees

IR populations

## [[1]]

## 
## [[2]]

## 
## [[3]]

## 
## [[4]]

## 
## [[5]]

## 
## [[6]]

Difference between scaled and unscaled IR populations

## [[1]]

## 
## [[2]]

## 
## [[3]]

% Difference of scaled and unscaled IR populations

## [[1]]

## 
## [[2]]

## 
## [[3]]